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Abstract 



We derive, from first principle, the Fokker-Planck equation describing the 

in 

■ non-perturbative real-time evolution of a gauge field at high temperatures on 

■ the spatial scale {g^T)^^ and the time scale (g^T)~^, where g is the gauge 



coupling and T is the temperature. The knowledge of the effective dynamics 
on such spatial and time scales is crucial for the understanding and quanti- 



Oh! 

1^ ■ tative calculation of the baryon number violation rate at high temperatures 



in the electroweak theory. The Fokker-Planck equation, which describes the 
diffusive motion of soft components of the gauge field in the space of field 
configurations, is shown to be gauge-invariant and consistent with the known 
thermodynamics of soft gauge modes. We discuss the ways the Fokker-Planck 
equation can be made useful for numerical simulations of long-distance dy- 
namics of hot gauge theories. 
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I. INTRODUCTION 



Despite a relatively long history, the problem of computing the rate of baryon number 
violation in the symmetric phase of the electroweak model still does not have the ultimate 
solution. Conventionally assumed to be (parametrically) of order O(a^T^), where T is the 
temperature, this rate has been re-evaluated to be suppressed further by an additional factor 
of a-w due to damping phenomena in hot plasma, which bring the rate down to 0(q;^T^) 
||T|J^. While the numerical coefficient of a^T^ is still not found, its knowledge is crucial for 
computing the baryon number produced by electroweak baryogenesis. 

The physics underlying B violation is sensitive to a particular type of modes in hot 
plasma: the almost static, soft magnetic modes with spatial momentum of order (y'^T.Q 
The same modes are well known to be the ones that cause the breakdown of perturbation 
theory at high temperatures. This flaw of the perturbation theory for hot gauge plasma 
has been viewed mostly from the context of static (equilibrium) thermal field theory, where 
the perturbative expansion of any equilibrium quantity, like the free energy, breaks down at 
the order that receives substantial contribution from static, soft g^T magnetic modes. The 
process of B violation at high temperatures presents another side of the complications due 
to soft magnetic modes: their real-time dynamics is also non-perturbative. 

For static characteristics, techniques have been develop to isolate and, in principle, com- 
pute the contributions from g'^T modes to physical quantities. In one of the approaches, 
the dimensional- reduction technique P], all hard, perturbative modes are integrated out 
analytically, and the theory for the remaining soft modes is essentially a 3d classical gauge 
theory at finite temperature. The problem of computing non-perturbative contributions to 
any physical quantity is reduced to solving the latter. 

In contrast, our knowledge of how to deal with the real-time dynamics of g^T modes in 



^Hereafter, by "soft" modes we will understand these modes. We assume the gauge coupling at 
the energy scale T to be arbitrarily small, g{T) <^ 1. 
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hot gauge plasma is much less reliable. It has been suggested that simulating soft modes in 
a classical theory can solve the problem in the quantum theory The justification for 

using a classical theory is the large occupation numbers of soft modes: when many bosons 
are on the same energy level, there is little difference between quantum and classical theories. 
Unfortunately, this method does not work because in a gauge theory, the dynamics of the 
soft modes depend crucially on the hard ones |l| . Even for classical theories, this dependence 
is present: the soft dynamics is different in classical theories with different lattice cutoffs p. 
The situation is in sharp contrast with the static case, where the hard sector only slightly 
(by an amount small in the weak-coupling limit) modifies the classical Hamiltonian of the 
soft one. 

At present time, it seems that the only way one can try to treat the real-time dynamics 
of soft modes is to "integrate out" the hard modes and see what is left out for the soft ones. 
The procedure should resemble Wilsonian approach to renormalization group. 

Integrating out hard modes in the real-time framework is not something new: the general 
formalism suitable for performing such integration, the "infiuence functional" approach, has 
been developed for quite a while 0J§]. Recently, the same method has been applied to the 
scalar 0^ theory In the latter case, the integration over hard modes leads to a stochastic 
equation for the soft sector, in which a damping term and a stochastic noise together keep 
the soft sector in thermal equilibrium. Technically, however, the effective equation is rather 
complicated: the noise and the damping are non-local in space and time, and both depend 
in a non-trivial way on the soft field. That the effective soft dynamics is complicated in such 
a simple model as the 0^ theory is one of the reasons that the same calculation has not been 
performed for the gauge field so far. 

One the other hand, to gain some insight into the effective dynamics of soft modes in 
hot gauge theories, the kinetic approach has been explored |T0|,3. In this approach, the 
soft sector is described a classical field, while the hard sector is replaced by a collection of 
classical particles moving on the soft background. Concerning the noise, it was suggested 
that it can be naturally introduced to the kinetic approach by considering the fiuctuations of 
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the distribution function of hard particles. The kinetic method yields a stochastic equation 
with damping and noise terms 0. It is remarkable that the damping and the noise are local 
in time (they are still non-local in space) and can be written is a relatively simple analytical 
form. The same is not true for the equation obtained in 0^ theory by the influence functional 
method. 

However, the use of the kinetic approach for hot gauge plasma is not well justified for our 
problem. While the kinetic equations are known to correctly reproduce the physical quanti- 
ties at gT momentum scale as the Debye screening length or the plasmon frequency ||TTl , p!2 



it is not clear whether these equations remain valid at the smaller scale of g^T. Moreover, in 
the framework of this approach, the stochastic noise remains a phenomenological quantity 
introduced by hand, and it is not obvious that the noise is Gaussian as assumed in Ref. 
or it has a more complicated distribution. Moreover, it has not been noticed in Ref. 
that there is an ambiguity in understanding the Langevin equation. Therefore, to reliably 
obtain the effective dynamics of the soft modes, one has to re-derive the equations from first 
principle, using a systematic approach. The influence functional technique seems to be the 
only available candidate for the latter. 

In this paper we present a systematic derivation of the soft dynamics in hot gauge 
theories, that is based the influence functional method. In this method, we divide the 
whole system in to a soft and a hard sector, and integrate over all modes of the hard one. 
Systematically keeping track of the order of the coupling g, after long calculation we recover 
the same set of Langevin equations as obtained from the kinetic approach. 

We furthermore observe that these equations are ambiguous, the fact that has not been 
noticed in the original derivation of . By fixing this ambiguity, we derive the Fokker-Planck 
equation that describes the diffusion in the space of field configurations in a unique way. We 
show that the Fokker-Planck equation is gauge invariant, a crucial test for our ideology if 
one takes into account that the operation of dividing the system into hard and soft sectors 
does not respect gauge invariance. We also point out that the Fokker-Planck equation is 
consistent with the known thermodynamics of soft modes, thus making a connection between 



the dynamics of the soft modes and their statics (or thermodynamics). The Fokker-Planck 
equation (Eq. (|i9|)) is the main result of our paper. 

This paper is organized as follows. In Sec. we briefly review the general-purpose tech- 
niques that will be used in the rest of the paper, namely, the Schwinger-Keldysh formalism, 
which provides the general framework for real-time non-equilibrium field theory, and the in- 
fluence functional method which allows one to derive the effective dynamics of a subsystem 
in a larger system. In Sec. |TT| we apply these techniques to the case of a hot pure gauge 
theory, and find the effective dynamics of the soft modes in this theory in the form of a 
Langevin equation. We discuss the ambiguity of the latter, and argue the way to fix this 
ambiguity. We then find the Fokker-Planck equation in the main text) that describes the 
diffusion of the gauge field in the space of field configurations, and check that, in the static 
limit, it leads to the same equihbrium Gibbs distribution of the soft modes as expected 
from the static method of dimensional reduction. We also check the gauge invariance of 
the Fokker-Planck equation. Finally, we present the concluding remarks in Sec. fV|, where 
the relevancy of the Fokker-Planck equation for numerical simulations of the plasma is also 
discussed. The Appendices contain various technical details. 

II. BRIEF REVIEW OF SCHWINGER-KELDYSH FORMALISM AND 
INFLUENCE FUNCTIONAL TECHNIQUE 

A. Schwinger-Keldysh (Close-Time-Path) formalism 

The most general framework for dealing with field theories in a real-time, non-equilibrium 
setting is the Schwinger-Keldysh, or close-time-path (OTP) formalism. In this section we 
will briefly describe this formalism, mostly to introduce the notations that will be used in 
the rest of the paper. Detail treatments of the Schwinger-Keldysh formalism can be found 



in Ref. 



Consider a quantum field theory of a generic field with some action S[(j)]. In the 
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Schwinger-Keldysh approach, the theory is, formally, put on a time contour which runs 
from some initial time moment, say, t = 0, to t = +00 and then goes back to t = (Fig. 
(|^)). The reason for having such a contour instead of the simple one running from t = —00 
to t = +CXD is the need of computing, in the real-time picture, the in-in matrix elements, 
i.e. those of a quantum operator sandwiched between two in-states instead of, as in the 
standard treatment of scattering processes, the in-out matrix elements between an in- and 
an out-state. The field in the upper and the lower parts of the contours need not to be 
the same, so we will denote in the upper and lower parts of the contour as 0i and 02, 
respectively, while reserving the notation for the field on the whole contour. The basic 
object in the Schwinger-Keldysh formalism is the generating functional, 

Z[Ji, J2] = J P02p[0i(t = 0), 02(t = 0)] exp(^i5[0i] - iS[(f)2] +ijdx (ji(x)0i(x) - 

-J2ix)Mx))) = Jv<j)p[<P{t = 0)] exp(zS[0]+a-0) (1) 
c 

where p[0i,02] is the density matrix of the system at t = 0, p[0i,02] = (02|p|0i) (we will 
take p to be that of thermal equilibrium), S'[0] = S'[0i] — 5'[02], J ■ (j) = Jcdx J(x)0(a;), and 
the index C of the last integral in Eq. (|Tp means that the variables in the integrations are 
taken on the contour C. Further we will largely ignore the boundary term p for simplicity of 
notations when it is not essential. All real-time physical quantities can be extracted from the 
generating functional Z; for example, the two-point correlation functions can be obtained 
by differentiating Z two times over J at J = 0, 
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c 



Z,' /l)00(x)0(y)exp(zS[0]) (2) 



where Zq = /^D0e*'^t'^l By putting x and y in different ways on the upper or the lower 
parts of the contour, the generic notation G corresponds to 4 different Green functions, Gu, 
G12, G21, G22, which in the operator language can be written as follows 

Gu{x,y) = {T(j){x)(j){y)), Gi2{x,y) = {(f){y)(j){x)), 

G22ix,y) = {f(f){x)(j){y)), G2iix,y) = {(f){x)(f){y)), (3) 
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where now is the field operator, T and T denote time ordering and anti-ordering, re- 
spectively, and (■ ■ ■) denotes the average over the thermal ensemble. Therefore, Feynman, 
anti-Feynman, and Wightman propagators can be recovered from Z. This list exhausts all 
interesting two-point real-time correlation functions. 

The Schwinger-Keldysh formalism also includes a set of Feynman rules, from which one 
can find any Green function in a weakly-coupled theory by computing the relevant Feynman 
diagrams. We will not go into details here, see [|l3l . 



For illustration and further reference, let us write down the explicit form of the thermal 
Green functions at equilibrium in some sample theories. The first example we will consider 
is the free massless scalar field. To compute the Green functions, it is convenient to make 
use of their operator form, Eq. (|^). Expanding the field operator in Fourier components, 

<P{t,x) = J (ape-IPl*+^P'' + ate^lPl^-P'^) 

where p = |p|, and noticing that, in thermal equilibrium, (aj,ap/) = (27r)'^2|p|(5(p — p')np, 
where Up is the Bose-Einstein distribution function, Up = (e'Pl/^ — 1)^^, we find the thermal 
Green functions for the free massless scalar case in momentum representation, 

. ( ^ + np Up \ 



1^2 _ p2 _|_ _ p2 _ 



np 1 + % 

,2_-r^2 I ,V ,.,2_„2 



cij^ — + te C(j"^ — p^ — te 



IpI 

G-l^{u,p) = ^^[il + np)5{u-\p\)+np5{u+\p\)) (4) 

(the superscript s means "scalar"). Let us notice a feature of the Green functions that 
will turn out to be useful in further discussion. When |p| is much smaller than T, one has 
Up ^ 1, and all four Green functions are approximately equal. G^^ ~ G^^ ~ G^^ ~ G^"^. 
This fact can be rephrased into the following statement: for soft modes (those with |p| <^ T), 



the path integral is dominated by such field configurations where the difference between 
the soft components of the field in the upper and lower parts of the contour is negligible. 

As the second example, we will write down the photon propagator for QED. The most 
convenient gauge to work with, in thermal field theory, is the Coulomb gauge where diA^ = 0. 
The photon propagator contains, beside the two contour indices a and P that can be either 
1 or 2, also two Lorentz indices. The spatial components of the photon propagator is 
proportional to the scalar propagator and has the same (spatial) Lorentz structure as that 
of the transverse projection operator, 

G%{u;,p) = (^6^^ - G:^(o;,p), = 1,2 (5) 

The 00 component of the photon propagator is the same for all values of the contour indices 
and does not have any singularity with respect to u, 

G°°(a;,p) = ^ a,/5 = l,2 (6) 

which refiects the non-propagating nature of the field component Aq. All other components 
of G vanish in the Coulomb gauge. 

In a non-Abelian gauge theory with weak coupling, the propagator of the gauge field, to 
the leading order of the coupling constant, is given essentially by the same equations as in 
the QED case, Eqs. The propagator now has two additional color indices a and b, but 

the dependence over these indices is trivial, i.e. via the Kronecker symbol 6"'^. 

B. The influence functional formalism and the stochastic equation 

In order to find the dynamics of soft modes in gauge plasma, we will divide the whole 
system into two subsystems, from which the first one contains soft modes with spatial 
momentum of order g'^T, while the second consists of hard modes. We then integrate over 
all the degrees of freedom in the second subsystem to obtain the effective theory for the 
soft modes. Before discussing the specific case of the gauge theory, we need to review the 
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general framework that allows one to find the effective theory describing a subsystem of a 
larger system in the context of real-time field theory. Integration over hard modes is usually 
performed in the Euclidean spacetime as, for example, in dimensional reduction approach 
0]. Doing the same in real-time, as we will see, is a little more complicated. The technique 
we will describe here is essentially the one developed and used in Ref. [§,0, where further 
details can be found. 

Let us consider some field theory which contains two set of degrees of freedom that will 
be generically denoted as and x- Eventually, and x will be identified with the soft and 
the hard modes in the hot gauge plasma, but for the moment we will regards them as just 
two subsystems making a full quantum system. The action of the full system is a sum of 
those for the two subsystems and an interaction term. 

We will try to find the effective theory of the field. More specifically, we will be interested 
in the n-point Green functions of the field 

G{xi...Xn) = Zq ^ JV(j)Vx (pixi) .. .(pixn) exp(iS[(f)] + iS[x] + iSint[<l),x]) 
c 

where 

Zo = J T>(f)Vx exp(iS'[0] + iS[x] + i5'int[0, x]) 
c 

and will derive a theory for alone, in which these Green functions are the same as in the 
original full theory. 

Take the integration over one finds 

G{xi ...Xn) = J T)(f)(f){xi) . . . (f){xn) exp(iS'[0] + iS'iF[0]) (7) 
c 

where we have introduced the so-called influence functional Sw, 

e^^'^t"^' =Jvx exp{iS[x] + ^5int[0, x]) (8) 
c 
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or, in component notation, 

^ J jy^^ exp(z5[xi] - iS[x2] + ^Sint[h, ^2, Xi, X2] 

So far all our equations are exact. Now we will make our first approximation by assuming 
that the integral in Eq. (|^ is saturated by such configuration of (f) where its values in the 
upper and lower parts of the contour are close to each other. 

Ad) Ad) 

01=0+^, 02 = 0-^, (9) 

where 

A0 < (10) 

In general, the condition (|1^) is not always satisfied. We make this assumption, keeping in 
mind the particular case of gauge theories, where is the soft modes and x represent hard 
ones. When the formalism is applied to the hot gauge theories, we will verify explicitly that 
(p!0|) is correct. However, not to claim any rigorousness, this assumption can be understood 
in the following simple way. Let us recall the field will be identified with the soft modes in 
plasma. We have seen (see Sec. that for the free scalar field the path integral is saturated 
by field configurations having the soft components on the two parts of the contour close to 
each other. Eq. (|l^) means that the same is valid for the soft components of the gauge field 
in the interacting theory. Again, this is not a rigorous proof and one needs to verify Eq. 
(p!0|) specifically for the case of hot gauge theory. 

Though 01 ~ 02, one cannot simply put 0i = 02 to the exponent in Eq. (^, since the 
exponential vanishes if 0i = 02. In fact, one can see that when A0 = 0, or 0i = 02, 5'[0] 
vanishes: the contributions to S[(j)] from the upper part and the lower part of the contour 
cancel each other. From unitarity it can be shown that S'if[0] also vanishes when A0 = 0. 
So, one should expand the exponent over A0, keeping the latter small but finite. For reasons 
that will be clear in subsequent discussions, we will be interested in terms of order A0 and 
(A0)^, but not (A0)^ and higher. With the symmetric definition of as in Eq. (|^), the 
expansion of S'[0], up to the (A0)^ term inclusively, is 
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S[cP] = S[<P^] - S[H = jdx ^A0(^) + O((A0)3) (11) 

The most general expression one may obtain expanding S'if[0] over A0, up to the order 
(A0)2, is 

Sw[4>] = - ldxj{x)A(P{x) + ^ldxdyA(l){x)N{x,y)A(P{y) 



(the time integration here, as well as in Eq. ([TlD , is along the Minkowskian time axis but 
not the time contour C) where we have introduced the following notations 

and 

t ( 5'S,F 6'S,p 6'S,p S'Sip \ 

^""'^^ 4 [6Mx)SMy) SMx)SMy) SMx)SMy) SMx)SMy) J 

Our conventions are chosen in such a way that in the case of gauge fields j and N are both 
real. Note that j and are, generally speaking, functionals depending on (p. 
Now the correlation function can be written in the following form, 

G(xi . . . x„) = Zq-i j V(t) VA(j) 0(xi) . . . 0(a;„) x 

X exp jdx [-^^ - 3{x)^ A(/)(a;) - ^ J dx dy A(j){x)N{x,y)A(j){y)^ (14) 

The integration over is Gaussian and can be easily taken. The result reads. 



G{x, . . . X J = Z^' J V<p 0(xi) . . . 0(x„) exp ("^^ " j) " ■ (^^ - ) 
The normalization factor Zq, analogously, can be written in the form. 



(15) 



Zo = JV(j) exp ^ 



U^s (16) 



50 

Eqs. (|1^,0) are the our final result and can be interpreted as follows. The Green function 
G{xi . . . Xn) of the full quantum system can be computed by by taking average over an 
ensemble of classical field configurations. Each field configuration 0(x) in the statistical 
ensemble has the following weight, 
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pW.),^Vexp(-lg-.).iV-.g-.)) 

Therefore, we have reduced the problem of computing the Green functions in the full quan- 
tum theory to the problem of averaging over a classical statistical ensemble with a known 
distribution of field configurations. In literature Eqs. (^,|I^) are usually interpreted 

in term of a stochastic equation. In this language, the dynamics of the field (j) is described 
by a Langevin equation with a noise term, 

where the stochastic source ^ is Gaussian distributed and is characterized by the following 
correlation function, 

{axny)) = N{x,y) (18) 

One should note, however, that N{x, y), in general, depends on cj), and the Langevin equation 
where the noise kernel depends on the solution is not very well-defined: to generate the noise, 
one need to know the solution beforehand, while the latter cannot be known without having 
generated the noise. The proper understanding of Eqs. ( ]T7|JT8| ) should be always Eqs. (|I5],|T6|). 



In further discussion, we will use the stochastic language of Eqs. ( p!7| , |l8|) , implicitly assuming 



that their meaning is given by Eqs. ([15| , [T6|) . 

The physical interpretation of j and can be seen from Eqs. (|T7| , [l8|) . First, according to 
Eq. (|18D, characterizes the amplitude of the noise coming from the x degrees of freedom 
that acts on 0. To counteract this noise, j must contain a damping that dissipates the 
energy coming from the noise. We will hence call j and the damping term and the noise 
kernel, respectively. 

The Langevin dynamics found above will serve as the starting point of our discussion 
of the hot gauge theory case. Let us recall here that the condition A0 ^ is the basic 
assumption in our derivation of Eqs. ( pTSD and (0), so one should explicitly verify that this 
condition is satisfied in the case of hot gauge theories. 
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III. SOFT MODES IN HOT GAUGE THEORIES 



Let us apply this formalism developed in the previous section to the case hot gauge 
theories. The Lagrangian of the theory is, 

where, beside the component notation A^, we use the standard matrix notation, 

= -tgAp^, F,, = -igF^^X 
F^u = d^A, + d,A^ + [A^, A,] (19) 

where are the generators of the gauge group G in fundamental representation. 

Consider the theory at some temperature T. We will divide the system into two, a "soft" 
and a "hard" , subsystems, and we want the soft subsystem to contain all non-perturbative 
modes. In Refs. 0,0] we have found that these modes have spatial momentum of order g'^T 
and frequency of order g^T. We also want most important hard modes, those with spatial 
momentum and frequency of order T, to be in the hard subsystem. Therefore, we take an 
intermediate momentum scale po = Ci9^T, and an intermediate frequency scale ujq = C2g^T, 
where ci and C2 are two arbitrary large numbers ci, C2 ^ 1 (which are still parametrically 
of order 1), and decompose the field into two parts, 

A^^ A^ + (20) 

where the new is the soft component which contains only modes with spatial momentum 
p < Po and frequency u < Uo, and the hard contains other modes. Note that the 
decomposition (^) does not respect gauge invariance, therefore one has to check the gauge 
invariance of the effective soft dynamics after it is derived. 
The field tensor is decomposed as 
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where the new F is related to the new A as in (|^), f^u = -D^Oy — D,^a^, where D^a^, is the 
covariant derivative of the hard field a on the background A, 

in which A^ is in the adjoint representation with the matrix elements A^f = —gf"'^'^A'^^. In 
general, we will use calligraphic letters for matrices in adjoint representation. 

To derive correctly the soft dynamics, it is crucial to know the typical magnitude of 
the field A^ and the soft field tensor F^^. These quantity can been estimated by using 
the fiuctuation-dissipation theorem We quote here the only the result. By order 

of magnitude, spatial components A1 ~ gT (or, in matrix notations, Ai ~ g'^T). This 
guarantees A to be non-perturbative: in the expression for the field tensor, Fij = diAj — 
djAi + [y4j,ylj], the last quadratic term is as important as the first two linear terms. The 
component Aq is usually set to by gauge fixing, however we will relax this condition by only 
demanding Aq to be of order g'^T, i.e. by a factor of g^ smaller than the spatial components. 
This allows one to perform the gauge transformation on A, A^ UA^U^^ + UdfJJ^^ ^ with 
U varying on the same spatial and time scales as A, i.e. (g'^T)"^ and {g'^T)~^ respectively, 
without breaking the conditions Ai ~ g'^T, Aq ~ g^T. The effective theory for A^ should be 
invariant under these "soft" gauge transformations. 

For the field tensor, the magnetic and electric components have different orders of mag- 
nitude. The magnetic components are F^j ~ diAj + ~ g'^T'^, while the electric ones are: 
Foi = doAi — diAo + [Ao, Ai] ~ g^T"^, where we have made use of the fact that (9o ~ ~ g'^T. 

Now let us apply the influence functional technique described in Sec. ||, identifying the 
soft fleld Afj, with (p and the hard fleld with x- To integrate over a^, one divides the 
action into the soft and the hard parts, 

S = S,[A,] + S^[a^,A^] 

where Sg contains only A, 

Ss[A,] = ^^iijdxFl, 
14 



while Sh contains the action for a and the interaction between A and a, 

S^a^, A^] = ^tr J dx (/^^ + 2F^4a^ a"] + ■ ■ ■) 

where dots denotes terms cubic or of higher orders on a. Following the general formalism, 
one should compute the influence functional, 

o 

To be able to integrate over a^, we have to fix the gauge for the hard field first. At finite tem- 
peratures, it is most convenient to work in the background Coulomb gauge, where a satisfies 
the constraint -D^a* = 0, since in this gauge only physical degrees of freedom (transverse 
gauge bosons) are propagating. One writes 

e5:p[A,] _ I exp(z5h[a^, A^] + logdet(D2 + [a\ Di])5iVia') (21) 
c 

where the commutator [a*, Di] is taken with respect to group indices. The determinant is 
usually rewritten into the form of an integration over the Faddeev- Popov ghost, but there 
is no need to introduce the latter in our context. 

Let us first show that logdet(D^ + [a*, Di]) is irrelevant. For this end we note that since 
the operator staying in the determinant does not contain any time derivative, it is equal to 
the integration of the 3d Euclidean log det over the contour, 

logdet(D2 + [a\Di]) = J logdet=^'^(D2 + [a\ D,]) 

c 

The 3d Euclidean determinant can be computed by calculating the one-loop Feynman dia- 
grams in the 3d theory, which is a quite complicated task. Fortunately, for our purpose we 
need only an order-of-magnitude estimation, which is 

logdet(D2 + a^A) ~ /rfx {g'F^^F^^ + g^f^^ 

The second term, g'^fijfij, is small compared to S'h[a,A] in Eq. (^Tj) which contains fijfij- 
The first term enters S'if[A] but is suppressed compared to S's[yl]. Therefore, the determinant 
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can be safely ignored. This result can be interpreted in the following way: the determinant is 
the ghost loop, but in the Coulomb gauge the ghost is not a propagating degree of freedom, 
therefore its contribution is the same as at zero temperature, which is suppressed by a power 
of the coupling constant. 

The influence functional, thus, can be written as 

e^^i^t^^l = J Va^ e^^'^['^^'^'*l(5(Aa*) = J Va^ VX exp (iS^a^, A^] + X'Via"^ (22) 

where we introduced an auxiliary variable A". According to Eqs. ([121) and ([13|) , to compute 
and N^y one need to take first and second derivatives of Sip over A^. 



A. Computation of the damping term j 

Let us find = 5S[A]/SA1^. One writes. 



r 



"2^/^''^^^ {W^jlf) (•5hk,^M]+^A'^^^««iexp(^^h[a^,A^]+^A'^V) 

= ^ /l^a^ VX (d^K, aT + K, FT + M"S^u) exp(z5h[a^, A^] + iX'^V^a'"') 

(23) 

where, for the simplicity of notations, we denote the sum over the upper and lower parts of 
the contour by the index "1+2". We will be interested primarily in the spatial components 
of . Running ahead, it turns out that the contributions from the upper and the lower 
parts of the contour are equal to each other, so one needs to compute only one of them and 
then multiply the result by a factor of 2. Let us concentrate on the contribution from the 
upper part of the contour (corresponding to the index 1). Using the constraint -DjO* = 
imposed by the integration over A, one can rewrite Eq. (PB|) into the form. 



/ = Zq^ j Va^VX {2D^[a^', a'] - [Doa°, A + [a°, - [a^ Aa^'] + [a\ A]) x 

X exp (iS^ [a. A] + iA" Aa"*) (24) 
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Let us introduce the propagators on the soft background, 

c 

HUx,y) = {a''\x)\\y)) = Z,' fva^a''\x)X\y) exp (zS^[a, A] + iX'^D^a'^') 



c 

(for avoiding cumbersome notations, we use lower group indices for G and H). In subsequent 
equations only the 11 component of these propagators will enter, so for the simplicity we 
will not write the index "11" explicitly. From Eq. (p^), one sees that j* can be expressed 
via these propagators as 



r{x)=gf 



abc 



2DlGZ{x, x) + Hl{x, x) + {D^,Gii{x, y) - D^,Gf,{x, y) - DlGt{x, y) 



(25) 



where we define the covariant derivatives of G{x, y) with respect to x and y as follows, 

d 

DlGab{x,y) = -^Gab{x,y) + A'"'{x)Gcb{x,y) 

DlGab{x,y) = ^Gab{x,y) - Gacix,y)A'\y) 

In order to find j, thus, one need to find the propagators G and H on the background field 
A. Let us specify the accuracy one needs to know these propagators. Since the Langevin 
equation has the form 55*8 [A] /(JA'^ = + we expect j to be of the same order as 
6S/6A'^ ~ D^A"^ ~ g^, therefore, one needs to find G and H up to contributions of order g"^. 

In Appendix ^ we showed that the temperature-dependent parts of (7°°, G^\ and H 
are of order g^ and thus can be ignored in Eq. (p^). Also, while G{x,x) is of order 1, 
its a;-dependent part is of proportional to g^ and varies on the spatial scale {g'^T)~^, so 
DiG{x,x) ~ g^. Therefore, Eq. ([25| ) can be simplified to 

r{x)=gr'^DfGii{x,y)\y^. (26) 

Instead of x and y, it is convenient to introduce the new coordinates X and s, 
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For simplicity, we will denote the propagator in these coordinates as G{X, s). As a function 
of s, G varies on the scale of T~^, while as a function of X the variation is over the spatial 
scale {g'^T)~^ and time scale {g^T)~^. In the coordinates X and s, Eq. ( p6D reads, 

= 9r'-(^ + A(A'))Gg(.Y,s) (27) 

To derive the last expression, we again have made use of the fact that dxG ~ when s —>■ 
{x — > y) and can be neglected. Thus, to compute j, one need to find the spatial components 
of G up to the order g"^ and plug into Eq. (|27|). 

In Appendix |^ the propagator G{X, s) is computed with the accuracy we need, i.e. up 
to terms of order g'^ inclusively. The result is somehow complicated, but some parts of it 
allow intuitive physical interpretation. First, there is a contribution to G that simply comes 
from gauge invariance. To extract this contribution from G, one writes, 

GZix, y) = U-\x, y)G'^{x, y) + G%{x, y) (28) 

where f/"^(x, y) are the matrix elements of the Wilson line connecting two point x and y in 
the adjoint representation. 



U''\x,y) = Te^^i-jdz^A^{z) 



y 



Go is the propagator in the absence of the background field (i.e., when A = 0), and G is 
of order g"^. All corrections of order g^, thus, is due to the Wilson line U"''^{x,y), while at 
the same time there is a part of the g"^ corrections that cannot be attributed to the latter. 
In particular, when x = y, or s = 0, propagator differs from Go (which is X-independent) 
by an amount of order g^. Substituting Eq. ( pS] ) into Eq. one sees that only the G 

term makes contribution to j. Since G is already of order g'^, to the order of g that we are 
interested in one can write, 



r(X) = gr"^-^Gli{X,s) 



It is convenient to perform the Wigner transformation on G, 

G{X,p) = Jdse'P'G{X,s) 

so j can be written in the form, 

r{X) = 2zgr^J^^p^Gli{X,p) 
Now let us quote the result of Appendix |^ for G (more precisely, on its 11 component) 



where 6n'^{X), with X = (t,x), is defined as follows, 

Sn^it, x) = -QjTT^ / du W^t, x; t - M, X - mv)v ■ E\t - m, x - uv) (30) 

and Gs is defined in Eqs. (H) and depends on p only via its module |p|. 

First let us substitute Eq. ( p9|) into the formula for j*. One notices immediately that 
for computing G^^ only the first term in the RHS of Eq. (|29| ) matters: other terms are odd 
under the permutation i ^ j. One finds, 

nX) = 2C,{G)g J ^ vW^iX) (31) 

where G2{G) is defined by f^<^dfbcd = C2{G)6''^ (for the SU(iV) group, C2(SU(A^)) = A^). 
Eq. (0) allows us to interpret 6n as the deviation of the distribution function from thermal 
equilibrium. Using Eq. (|30|), j can be rewritten via the soft field in the following way, 

fUt, x) = -2C2(G)g^ [ Jdu -^f/^Vt, x; t - M, X - UY)vVE^Ut - m, x - mv) 

J {2n)'^ J a|p| 

One can now write dp = p^ d\p\ dv, where v is the velocity that corresponds to p, v = p/|p|, 
and take the integration over d\p\. The result reads, 
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\2tx J 

Now in the integral it is convenient change the variables u and v to a new variable y = x— mv, 
and one finds, 

fit, X) = ^Mgg! hy ,r\t, x; t - |x - y|. y) - ,7 E^Ht - |x - y|. y) (32) 

IzTT J |x — yp 

It is natural to assume that the integral in Eq. (|30|) is saturated by such y that |x — y| ~ 
{g^T)~^. In fact, for larger distances between y and x the Wilson line becomes a rapidly 
varying function, so the contribution from theses distances vanishes. Since |x— y| ~ (g'^T)^^, 
one can replace, in Eq. (|32D, t — |x — y| by t, as the typical time scale of the evolution of 
soft modes is much larger than (g'^T)^^. Therefore, the dependence of on the soft field 
can be considered as local in time (but still nonlocal in space) , 

r(i,x) = ^M£l9!l! Ut'-(t,x;i,y)(i^i^):(i^£'i(t,y) 

izvr J |x ^ y\ 

One can substitute = —A^ — DiAo to this equation. The term with Aq can be taken by 
integration by part and the result vanishes due to the identity 

D: (^[/°^(t,x;t,y) ^^~^^^^^^~^^' ) = (33) 

so one obtains finally. 

Note that we have computed j by multiplying the "1" component in Eq. (|23| ) by 2. One 
can check that the "2" component, which can be found by the same method as described 
above, is equal to the "1" component. 

B. Computation of the noise kernel N 

The noise kernel, according to the general formula (0), is equal to 
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(35) 

11+12+21+22 



A6A1{x)6A]{y) 

where we have indicated that the sum over all possible locations of x and y (on upper or 
lower parts of the contour) should be taken. Similar to the case of j, it turns out that all 
four terms corresponding to x and y on different parts of the contour C are equal, so one 
need to compute only one of them and then multiply the result by 4. Let us pick the 12 
component and write, 

^^^"'^^ = -^ Mf.(x)M^,(,) ^"[^] 

In the following we will not write the indices 1 and 2, implicitly assume that x is on the 
upper and y is on the lower parts of the contour. Using Eq. ([2^), one finds. 



NZix.y) = {[a''ix),D,a'{x)r[a\y),D,a'iy)]') - {[a\x), D,a'{x)r) {[a\y), D,a\y)]') = 

= r'f'f {Gl{x, y)d'td]G',i{x, y) + 9|G5(x, y)d^Gl{x, y)) (36) 

where we have replaced the covariant derivative D by the simple derivative d. It is possible 
to make such replacement since we will be interested only in the leading order. In fact, 
D = d + A, and for hard field like a, d ^ T while A ~ g^T. 

Since we are interested in the case when x — y ~ {g'^T)~^, we need to know the propagator 
in the regime of such separation between x and y. To the leading order on g, the calculation 
of the propagator for x — y ~ {g'^T)~^ is carried out in Appendix The result, for the 12 
component, reads 

G'l^^^^y) = /(2^ {fMf,Hy)y\5- - i^)(npe-ipi(-*')+-(--) + 

+ (l + np)e*IPl(*-*')-*P(''-y) 

where x = (t,x.), y = (t', y), and /p(x) is, for each value of the coordinate x, a group element 
in adjoint representation which is defined as follows, 

f;\x) = U^\t,^;0,^-^^t) (37) 
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(as before, v = p/|p|) provided the background A is turned of for t < 0. Now let us 
substitute this propagator to Eq. (^). The first term in N is 

dpdp' ■ f ■■ p'^p'\ / • • p'^p'^\ 



p/2 



X 



X {U{x)f;\y)Y {fA^)f,My)y (npe-^IPl(*-*')+-(^--) + (1 + np)e^lPl(-*')-P(^'-)) x 

X (np'e-^IP'l(*-*')+*P'(^-y) + (1 + npOe^'P'l(*-*')-^P("~^^) (38) 

Notice that p, p' are of order T, while t — t' and x — y are of order {g'^T)~^, therefore the 
integrand in Eq. (^) is a rapidly varying function unless p and p' are close to each other. 
The integral, thus, is dominated by this region. Denoting p — p' = Ap, and notice that 
when Ap is small one has |p| — |p'| = v ■ Ap, one can write Eq. (^) into the following form, 

rj'''G'Je{x,y)d!d]G%{x,y)=4f^ X 

xnp(l + ^p)py e-^^-'^P^^-^'^+^^P^^-y) 

Now one can take the integration over Ap, which yields a delta function, 

dp I j: I \ /.— 1/ \\'^^ I £ I \ £—\l \\'^^ 



r^^f^fG':t{x,y)d:d]G',j{x,y) = 4^'^^ / (^^^(/p(^)/p ^(z/)) {Ui^)f^\y)) x 

xnp(l + np)p'p>6{x - y - v(t - t')) 

Due to the the delta function, the integrand is non-zero only when x — ?/ = (t — t', x — y) is 
proportional to f = (1, v). In this case, one can verify using Eq. ( |57| ) the following identity, 

fp{x)fp\y) = Uix,y) 

which is independent of p and can be taken out of the integration. Therefore, 

r'f'fGf,{x,y)d'[dy^G%{x,y) = f^fU^%x,y)U'^ {x,y) x 

X ^Vnp(l + np)5(x - y - v(t - t')) 

The pre-integral factor can be written as —iiT"'U{x,y)T^U~^{x,y) where T° are the gener- 
ators of the gauge group in adjoint representation, (T")'"^ = —f"-^^. To write this factor in 
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a simple way, one notice that the following identity holds for any matrix U in the adjoint 
representation of the group G, 

The second term in Eq. ( |5BD can be computed analogously, and the result is equal to that 
of the first term. One obtains, 

iV^^(t, x; t', y) = 2U^\x, y) J ^t;V%p(l + np)5(x - y - v(t - t')) 



Taking the integration over dp, one finds, 

^3 = ^J^C'C - <' - |x - y|) + ^« - *' + |x - y|)) 



127r |x — y| 

(39) 

Now since x — y ~ (g'^T)^^ which is much smaller than the time scale interested in, one can 



ignore the time non-locality in and replace Eq. (39) with a local one 



Nab[t,yi;t ,y) = U (t,x;t,y) ■ -j d{t-t) (40) 

DTT |x — yp 



Using the same method, it is easy to check that the 11, 21 and 22 components in Eq. ( |35D 
is equal to the 12 component, so can be computed by taking the 12 component and 
multiplying to 4, as we have done. Therefore, Eq. (BO) is the noise kernel we want to find. 



C. Verification of the initial assumptions A A <^ A 

To make sure that we have derived the right equation, one has to verify that the basic 
assumption used in the derivation, namely, that the difference between the values of the 
soft field A in the upper and lower parts of the contour is small, AA ^ A, is valid (more 
precisely, that the path integral is saturated by field configurations with A A <^ A). For this 
end, let us recall that the integral over AA is taken with the weight. 



exp (^-i I dx dyAA%x)Nab{x,y)AA\y) 
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So, knowing the order of magnitude of A^, one can estimate how large is a typical AA. As 
we have found before, N{x, y) is nonzero when |x — y| ~ a;o — yo ~ {q'^T)'^, and for these 
values of x and y the typical value of is, according to Eq. (|39D, 



N{x, y) ~ /T^ ■ ■ ~ g'T' (41) 

(x - y f [xq - yo) 



Now to find the typical value of AA, one writes 

dxdyN{AA''Y 



provided the integral is limited to a region of x and y with linear size of order {g'^T)~^. 
Given the typical value of A^ in Eq. (0), one obtains the estimate A A"- ~ g'^T. Recall that 
the typical value of the soft field is A"' ~ gT, one sees that AA <^ A. 

Now we will explain why we have treated the terms linear and quadratic on AA on equal 
footing while neglecting cubic and higher terms in the expansion of Sip. Naively, if AA <^ A, 
the linear term j ■ AA is much larger than the quadratic term AA ■ N ■ AA. However, we have 
seen that there is a delicate cancelation in j: the leading and next-to- leading orders gives no 
contributions to j, the first nonzero contribution comes from the order g"^. In contrast, there 
is no such cancelation in A^. More careful order-of-magnitude analysis show that these two 
terms are equally important. Since there is no cancelation in the second term A A ■ N ■ A A, 
higher terms are certainly more suppressed than the second-order one and can be neglected. 

D. Fixing the ambiguity of the Langevin equation: the Fokker-Planck equation 

Above we have computed j and A^, which are the two quantities one need to know in 
order to write down the Langevin equation for the soft modes. The Langevin equation has 
the following form. 



127r J |x — yp 



(r'(i, x)«''(i', y)) = (£_l):(£_l)!t,-(i, , _ ,) (43) 

OTT |x — yp 
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Notice that in the LHS of Eq. (H) should stays 5S[A]/6A^ = -Do^oi + DkFki, but Do^oi is 
smaller than D^F^i by a factor of g'^, so the former can be neglected. These equations have 
been first written down in Ref. 0. However, an important point has been missed there: the 
Langevin equation with local noise and local damping, where the noise kernel is a function 
of the coordinate is ill-defined, in the sense that different discretization of this Langevin 
equation leads to different evolution in the continuum limit. Therefore, one has to be more 
careful in interpreting Eqs. (^ , ^31) as the equation describing soft dynamics in hot gauge 



plasma. 

The ambiguity of the Langevin equation with local noise and damping depending on the 
coordinate is known in literature but since the topic is important for our discussion. 



let us illustrate it on a simple example of Brownian motion of a one dimensional particle 
describing by the equation, 

-f{x)x + V'{x) = C (44) 



Note there is no x term as in Eq. (|4^) , which means that we are in the regime of over- 
damping. The coefficient 7 represents friction and varies with x. The noise kernel will be 
chosen to be local in time, but also x-dependent, 

mm) = '^T^ix)5it-t') (45) 

The unique way to characterize the Brownian motion is via the Fokker-Planck equation, 
which describes the diffusion of the particle. Denoting p{t, x) to be distribution function 
of particle (the probability of finding the particle at the point x at time moment t), the 
Fokker-Planck equation is the one describing the evolution of p(t, x) with t. 

In Appendix |C| we demonstrate the ambiguity of the Fokker-Planck equation for this 
simple one-dimensional case. We check that for two different ways of discretizing Eqs. 
(0,^), one obtains two different Fokker-Planck equation. For example, if one approximates 
Eqs. ( PI^SD by the following finite-difference equations 
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and then takes the hmit At 0, then the Fokker-Planck equation is 

while, if one chooses to compute 7 at x{t + At) instead of x{t), i.e. use the following dis- 



cretization, 



/ / , , , x(t + At) — sft) ^, , 
7(x(t + At))^ ^ = iit) 



. ^^ 2r7(x(t + At)) , 

(e(t)e(t + nAt)) = — ^-^^ — ^5„o 



then the Fokker-Planck equation becomes 

dp ^ d f 1 d 

dt dx dx ^ ' ^ I ' dx \ 'J 



^ = 711^ 7 p)]+7^\—p] (47) 



It is clear that this equation is different from Eq. (46) which has been derived using another 
scheme of approximation. In general, by choosing different discretization scheme, one could 
end up with the equation 



dt dx\ dx^ W 9a; y 7 

with arbitrary value of a. Eq. ( ^Uj ) corresponds to a = 1 while Eq. (^7]) corresponds to 
a = —2. Which value of a one should take for the Fokker-Planck equation in the case of 
gauge theory? 

We argue that the value a = is the correct one. Note that both discretizations of Eqs. 



(0,^) are asymmetric: the damping and the source kernel are evaluated at either end of 
the time interval (t, t + At). It seems that the more "correct" choice would be the one where 
these quantities are evaluated at the middle of this interval. In Appendix we consider 
such a scheme. In this scheme, we add a small second derivative term ex to the Langevin 
equation and discretize it in the following way, 
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{((t)((t + nAt)) = 2T',(x(t))S„o 

The small second- derivative term must be added for the prescription to have a limit at 
At 0. Moreover, in the case of gauge theory the second time derivative of A does 
really present in the equation (we ignore it in Eq. (^) since it is much smaller than spatial 
derivatives). The limit At — > should be taken before the limit e ^ 0. In Appendix we 
show that the Fokker-Planck equation, for this prescription, has the form, 

dp(t,x) d ( T d , \ d (V'ix) , \ 
which can be written as 

^ + ^ = 
dt dx 

where j = — — is the probability current, which is a sum of the current coming from 
diffusion (the first term) and from the overall drift due to the force V (the second term). 
The static solution to Eq. (^Sf ) is given by the barometric formula, 

p(t, x) ~ e T 

One can check this by verifying that for this probability distribution, the probability current 
j vanishes. 



Eqs. (^,iD is the mult i- dimensional version of Eq. ( ^^5| ) with a particular choice of 
7. It is straightforward to generalize Eq. (^Sf ) and write the Fokker-Planck equation for our 
case. Let p{t, A) be the distribution function, then 

|p(,A) ^ ^x..-^ (h-)-(x.y; A) (r^^^^ . I^M'-A))) («) 

where 7"^ is the inverse of the operator 7, 



7t {^,y;A) = t/"''(x,y; A 



12tt |x — y 
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and V{A) = | jd^F^-F"'- is the static energy of the field configuration A. 

The Fokker-Planck equation, Eq. ( p9D , is the central result of this paper. 

While the Langevin equation (^,^) is ambiguous, the Fokker-Planck equation ( ^Qf ) 
defines the evolution of soft modes in a unique way. There is still a problem with the 
degeneracy of 7, which is related to the gauge invariance of the Fokker-Planck equation. We 
defer this question to the end of this section. 

Similarly with the one- dimensional case, Eq. (^91) possesses a static solution, 

p(A) ~ exp (-^) (50) 

Therefore, the Fokker-Planck equation we have derived is consistent with the thermody- 
namics of the soft modes: it is known from dimensional reduction that the statics of soft 
modes is described by a classical ensemble with the barometric distribution as in Eq. (|50|). 

We did not present the rigorous proof that Eq. (|^) is the correct choice of the Fokker- 
Planck equation. To rigorously show the latter, one need to go one step backward from 
Eq. (^,^), i.e. not neglect the non- locality in damping and noise kernel and carefully 



study the behavior of the solution on time scales much larger than that of the non-locality. 
Unfortunately, this program seems to be very complicated. We believe, however, that the 
arguments presented above convincingly showed that Eq. (|49|) is the real Fokker-Planck 
equation describing the soft modes. 

Let us finish this subsection by noticing the Markovian character in the evolution of the 
soft modes. The time evolution of A can be visualized as a large numbers of transitions 
t^t + /S.t^t + At2 — * ■ ■ ■ where each transition occurs independently from others. 
We emphasize that the Markovian character is manifest only at large time scales (such 
as (g'^T)^^): on the time scale of {g'^T)^^ the evolution of soft modes is certainly non- 
Markovian. However, there is no substantial evolution during such small time intervals. 
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E. Gauge invariance of effective soft dynamics and degeneracy of 7 



In this subsection, we will show that the Fokker-Planck equation we have found is gauge 
invariant. Recall that for A we did not fix the Aq = gauge, but only impose a rather 
relaxed condition that Aq ^ g^, which allows one to make any soft gauge transformation with 
parameter varying on the spatial scale {g^T)^^ and time scale of order (g'^T)"^. Therefore, 
the effective dynamics of the gauge field must not break the invariance under these gauge 
transformations. The verification of this fact is a crucial test for our formalism, since, as we 
already noted, the division of the system into soft and hard modes is not a gauge invariant 
operation and it is important to check that the final result is nevertheless gauge invariant. 

Let us first define what is is soft gauge invariance in term of the distribution function 
p(A). The gauge invariance of p simply means that p(A) is the same for equivalent field 
configurations A, 

p(A) = p(A^) (51) 

where Ai = UAiU~'^ + UdiU~^. We will now check that the RHS of the Fokker-Planck 
equation is also gauge-independent, so if p satisfies Eq. (|^) at some time moment it will 
satisfy at later times. 

Under gauge transformation, 5/5A transforms as 

5AU 5A5AU 5 A 5 A 

where U is the matrix U in the adjoint representation. The matrix 7(x, y; A) transforms as 

7(x,y;A^) =W-i(x)7(x,y;A)W(y) 

so 7~^ should transform analogously. One sees that the RHS of Eq. (^Uj) is invariant under 
gauge transformations. 

However, the real situation is more complex. In Eq. (^91) we implicitly assume 7 is not 
degenerate in writing 7^^. However, it is easy to see that 7 is degenerate! In fact, assume 
one wants to find /3j = 7jj^aj. For that one needs to solve the equation, 

29 



Jdy -fij (x, y ; A)/3j (y ) = (x) (52) 

Taking the derivative Di on this equation and summing over i, making use of the identity 
Df7jj(x, y) = (see Eq. (|33|)), one finds -Djaj = 0. Therefore Eq. ( |5^ ) possesses a solution 
only when Diai = 0. 

Moreover, provided there exists a (3 that satisfies Eq. (0), one can easily check that 
A' = A + DiT], where rj is an arbitrary function, is also a solution. Therefore, we conclude 
that 7 is degenerate: 7"^ can be defined only when it acts on a function satisfying the 
requirement D^ai = 0, and the result itself is not uniquely defined but only up to a full 
covariant derivative. 

Then how one can understand the Fokker-Planck equation (^)? One will see that this 
equation is still well-defined just because of the gauge invariance. 

Let us verify that one can act 7^"*^ on |^ + |jp. One has to check that 

It is easy to see that ~ DiDjFij = 0. Therefore, what really needs to be checked is 

Dij^ = 0. This condition comes from the gauge invariance of p. In fact, under infinitesimal 
transformations, Eq. ( pT| ) reads, 

p{Ai) = p{Ai + DiO) 

for any function a. This is equivalent to 

dyi . Aa(x) = 
dAj(x) 

Integrating by part this equation, and making use of the arbitrariness of a, one finds, 

which is exactly the condition for 7~^|f to exist. Therefore, we have checked that one can 
define 
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/f (X) ^ (t^I . i^p) (53) 

in the sense that 7/ = + |^p, if p satisfy the gauge invariance condition (pl|). 

There is still an ambiguity in the definition of /f (x): this quantity is defined only up to 
a full covariant derivative. One may ask whether different choice of this quantity lead to 
different evolution of p. We will see that this is not the case, provided one defines 7"^ in a 
gauge invariant way. The latter means the requirement that / transforms covariantly under 
gauge transformations, 

/(x;A^) = [/(x)/(x;A)f/-i(x) 

After imposing this condition, / is still not defined uniquely. One can replace 

/f (x; A) ^ /f (x; A) + D,c\^- A) (54) 

where the only requirement on c(x; A) is that it satisfies the equation 

c(x; A^) = f/(x)Ac(x; A)t/-^(x) (55) 

for any function ?7(x). After the replacement (^), the RHS of the Fokker-Planck equation 
(^) picks up a term proportional to 

For the Fokker-Planck equation to be uniquely defined, one needs to show that this expres- 
sion vanishes. 

Let us consider Eq. (|55D for infinitesimal gauge transformations U = 1+a, Ai Ai—Dia. 
It reads, 

Ac(x; A) - [A«(x), c(x; A)] - f dy ^f^)^^ D,a{y) = Ac(x, A) + [a(x), Ac(x; A)] 
Taking the integral over dy by part, one finds 



A^/^y^J^^«(y) = Z)n«(x),c(x; A)] 
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which imphes, 



"(y) = ["(x),c(x; A) 



In components, this equation reads 



dt 6A\{y) Mf(y) 
Since this equation must be vahd for any function a, one finds, 

d ScH^- A) ,gfc,^c^y^S^^ ^ -z^rV(x)5(x - y) 



Now putting a = b and take the sum over a, and taking y — x, we obtain 



d 5c"fx:A) 



dy' 6Af{y) 



+ zgr'Ali^ 



5c'^(x; A) 
' SAfi^) 







Taking the integral over x, and noticing that 



dt SAfiy) 
is a full derivative, one finds, 

'59iC'^(x; a: 



+ 



d (Jc'^fxiA) 



dx' 6A^{y) 



y^x 



dx 



igr''Al{x) 



5cHx, A) 







The LHS of this equation is the same as the expression we want to show to be vanished, Eq. 
(^), thus we have demonstrated that the RHS of the Fokker-Planck equation is independent 
of the way we choose 7"^^. So, despite the degeneracy of 7, the Fokker-Planck equation is 
still well-defined and gauge invariant. 



IV. CONCLUSION 

In this paper we have derived the Fokker-Planck equation describing the non-perturbative 
real-time dynamics of the soft modes in hot gauge plasma, Eq. (^). With some reservations, 
the same dynamics can be described by the Langevin equation, Eq. (^,|^). The dynamics 
is diffusive in nature, where the word "diffusive" refers to the evolution in the space of 
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field configurations: assuming a given configuration of the soft field at some time moment, 
the field will "diffuse" and eventually fill the whole space of field configurations with the 
thermal distribution during a time interval of order {g'^T)^^. The behavior is totally different 
from that predicted by the classical field equation, in contrast to the case of scalar theory 
T3| where there is a time regime where the classical field equation can be applied. The 



difference between the behavior of the gauge field and of a scalar field originates from the 
fact that in the gauge case the soft modes are in the regime of over-damping |l|] while in 
the scalar case they are not. We believe that the understanding the diffusive nature of 
the soft field evolution in gauge theories is crucial for any discussion of physical processes 
involving real-time dynamics of soft modes, for example, baryon number violation at high 
temperatures. 

While from Fokker-Planck equation is the fundamental equation of long-distance real- 
time dynamics at high temperatures, from the point of view of numerical simulations the 
Langevin equation (^,^) may have its advantage. However, even the latter is hardly 
suitable for real simulations due to the spatial non-locality. Moreover, the Langevin equation 
is ambiguous. At this moment, it seems that the only practical way is to simulate a classical 
theory whose soft modes have the same dynamics as the hot quantum theory. Some proposals 
have been put forward in this direction. One is to simply simulate a classical lattice theory of 



a gauge field |jT6[ , another is to simulate the kinetic of a plasma consisting of a soft field and 
a collection of hard particles [0. The latter approach might be more economical, though 
the rigorous proof that the dynamics of the soft modes in this picture can be made identical 
to that of the quantum theory is still lacking. The Fokker-Planck equation derived in this 
paper can serve as a criterion for choosing the right simulation procedure, which, hopefully, 
will leads to the ultimate solution to the problem of hot baryon number violation rate. 
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APPENDIX A: HARD PROPAGATORS ON SOFT BACKGROUND: SMALL 

SEPARATIONS 

In this Appendix, we will compute the propagators of the hard field a and the auxiliary 
field A on the soft background A, G (aa) and H ~ (aA). We will need to find these 
propagators up to the correction of order g'^ inclusively. Since A ~ g^T, this means that one 
need to find also correction of the second order with respect to A. 

We will be especially interested in spatial components of G, 



We will first compute this propagator at tree level, i.e. neglecting the non-linear (on a) terms 
in S'h, and then show that the result does not change when one turns on the interaction 
between hard modes. 



The convenient way to compute the propagators at tree level is to introduce a source 
j coupled linearly to a and study the linear response of the system in the presence of this 
source. In this way, the propagators can be computed using the following formulas. 




1. Tree level 
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where (• ■ ■)j denotes the mean value of in the presence of the source j, 

(a'^'^), = Z^^\]\ J Vaf'VX a'^'^(x) exp {iS^[a, A] + iX^'D.a"' - ij'^a"'' 
{X^)j = Z^^[j] J Va" VX A"(x) exp {iSh[a, A] + iX^Dia^' - ij^a^^) (Al) 

where 

Zo[j] = J Va^ VX exp (iS^la, A] + iX^D. 



At tree level, finding these main values is equivalent to finding the saddle point of the 
exponent in Eq. (^), i.e. the solution to the following set of equation, 



D^a" - Df^iD ■ a) + 27^" - b^'DiX - = 

Aa* = (A2) 



where D'^ = D^j^D^^. We will need to know the solution to Eqs. to the accuracy of 

(we will see that this accuracy is necessary for finding the propagators with the accuracy of 
(7^). First let us find aP . Substituting /i = to Eq. (|A^), one obtains. 



where denotes the field tensor in adjoint representation, T°'^ = —gf"'^^F'^, from which one 
can express a° via other variables, 

a' = -B~Y + 2B-'{W) (A3) 

Notice that a* is hard, and J-'oi ^ g^, so without making error of order 0{g^), this equation 
can be written as 

= -D-2j0 + 2J^oid-^a' (A4) 
Since G^^ ~ da^/dj^, and J^oi ~ 9^, one finds immediately 

G-oo = + 0((7^) = G°°|t=o + 0(/) 
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where we have made use of the fact that the leading term in G^" does not contain time 
derivative and thus is temperature-independent. From Eq. ( |A4| ) one can also estimate G^*. 
Indeed, ~ da^/dj^, and from Eq. ( |7^ ) a° depends on j* only through a*. Since Fqi ~ g^, 
one finds, 

Let us turn to the equation for a\ One writes, 

D^a' + Di{DQa^) - 2Ti,a' + AA = f 

Substituting aP found in Eq. (|A3|) , one obtains 

D'^a' + Di{-DoD-'^f + 2J^Qjdod-'^ai) - ^Ti^a? - AA = 3' (A5) 

Let us express A via the other unknowns. For this end we apply Di to Eq. (|A5| ) and make 
use of the constraint -Dia* = 0. One finds, 

- {V,T,,)o? - D^AD-^jO - D^A = Aj" (A6) 

Denoting l^iTij = J^j, which is of order g^, one can extract A from Eq. (|X6|), 

A = -DoB-'f - B-\Da') - Td-^d 

From this one also finds 

HUx,y) = HUx,y)\T=o + 0{g') 

i.e. the temperature-dependent part of A is of order g^. Substitute a° and A to the equation 
for a\ the latter reads, 

D^a, + 27^kdod-^d,ak - 2^,fca'= + Jkd-^d.a!' = f - D.-D-^D^f + 0(/) (A7) 

To find the Green function, now we adapt the technique for finding the kinetic equation (see 
to the case with background field. First, according to Eq. (|A7|), the spatial components 
of the Green function G satisfies the following equations 
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+ Ax{x)] Gij{x,y)+i2J^ok{x)^ + Jk{x)] —Gkj{x,y)- 



The same equation can be also written as 

-2G,kix,y)J'kj{y) = -iSljix,y) + 

where 5*j(a;, = 6ij6{x,y) — DiD^'^Dj6{x,y). For simplicity wc did not write the contom 
indices (1 and 2) in our equations, with the imphcit understanding that one can put these 
indices back in at the end of the calculations. Let us subtract the two equations from each 
others, and divide the result by 2. The RHS vanishes while there are 3 different contributions 
in the LHS which we will denote as 

(ffj (^^o.(y) + Imv?! (A9) 

h = -{Fik{x)Gkj{x,y) - Gik{x,y)Fkj{y)) (AlO) 
The obtained equation reads, 

/i + /2 + /a = (All) 

First, let us consider /i. It is convenient to use, instead of x and y, the new coordinates 
X and s which are related to the old coordinates as 
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For simplicity we will denote the propagator G{x,y) in the new coordinates as G{X, s). As 
a function of s, G{X, s) varies on the scale of T~^, but G is a slowly varying function of X 
which changes on the spatial scale of (g'^T)^^ and time scale of ((^^T)"^. We will evaluate 
each contribution Ji, I2 and J3 to the accuracy of 0{g^). First, consider Ji. Expanding Eq. 
3), making use of the fact that A ~ g^T, dx ~ g^T, ds ~ T, one finds, 



h 



d dC 



+ 



A 



dG 

ds,, 



If dG 

+ 2 r^'sx: 



l{s^d.A,,§-] + l{s^A,G} + 



+- 



1 



A^G +-{A^,s''d,A^}G 



(A12) 



where we [■ ■ ■] denotes commutator, {■ ■ ■} denotes anti-commutator, and all A's are computed 
at the point X (from now, when A and F are written without arguments we will implicitly 
assume that they are computed at the point X). 

To the zeroth order on g, the propagator is equal to that in the absence of the background 
field. Gab = GoS""^, where Go is defined in Eqs. (j^) and (|]). When the background is 
present, Go{x, y) is not a gauge invariant quantity. To correct this, one can include a Wilson 
line U{x,y) to GQ{x,y), and the result will be gauge invariant. However, there may be 



contributions to G which cannot be attributed to the Wilson line, so we write. 



Gab{x, y) = Go{x, y)U'^ (x, y) + Gah{x, y) 



(A13) 



Expanding the first term in the RHS of Eq. ( |A13D on g^, one obtains 



G = Go - (sM^)Go + ^(sM^)2Go + G + 0(^« 



(A14) 



So far, we did not make any assumption on G. All we know about G is that it is a correction 
to Go, thus it has order of magnitude of g"^ or less. We will assume, and verify subsequently, 
that G in reality has the order of magnitude of g'^. Therefore, the g^ correction in G comes 
entirely from the Wilson line, while the same is not true for the g'^ correction. 

Substituting Eq. (|A14|) to Eq. (|A12|) and neglecting any contributions of order or 
higher, one finds, 

^ dG ,, ^ dGn 1 ,, r ^ \ < dGn 
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where = djdX^ + [^^(X), . . .]. Now will will split l\ into two parts and evaluate each 
separately. The first part is 

dG 9Gn 



while the second one is 

ht = -s''{J',.,s'A,}^ (A15) 



It is useful to rewrite Iia in the Wigner representation, i.e. in the momentum representation 
with respect to s. Introducing the notation, 

G{X,p) = Jd\se'P'G{X,s) 

and similar formulas for other Go and G, one finds 

ha{X,p) = -rp^^D,G{X,p)+p^T,.{X)^^^ 

OPu 

Now, let us take into account the fact that the structure of Go with respect to group and 
Lorentz indices is 

{G,)Z{P) = {s^' - ^) (A16) 

where Gs{p) is the propagator of the free massless scalar field. After some calculations, Iia 
can be expressed via Gs as follows. 



ha = ( S'' - ^] P^-^/^-^ + ^i^o^p' + J'o.vlGs + ^(^fc^F'' + J'kjP')Gs (Al7) 



Now let us turn to In,. From Eq. (|A15|) one sees that to compute In, to the order of 



(yf^, one need to consider only the case when /x, v and A are all spatial indices. It is more 
convenient to leave In, in the coordinate representation, where 

One can substitute Go to this equation. For our purpose, it is sufficient to notice that in 
the coordinate representation, the most general tensor form of Go is 
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where and are functions depending only on the module of s. In term of these functions, 
lib can be rewritten as 



^sV {{F,,,Ai}s^ + {F^,,Ai}s') G2 



So, we have computed Ii. For I2, substituting Eq. (|A16|) to Eq. (^), we find, in the 
momentum representation, 

h = -^(W + :Fo,p')Gs - T^ip'J, - TpJi)Gs (A18) 



Notice that the first term in the RHS of Eq. ( |A18| ) cancel the second term in Eq. ( |A17| ). 



Now consider I^. Making use of Eq. (|A14|) , can be divided into two parts. Is = /3a + -^36 
where 



'3a 



■[^,Go] 



\[s^dk:F,G^] + [:F,{s^Ak)G 



where the matrix multiplication is understood as occurring for both (spatial) Lorentz and 
group indices and both A and JF are computed at the point X. I^a is of order (7^, while I^b 
contains 0{g^) contributions. For I^a, the calculations is quite simple. Using Eq. (|A16| ), one 
obtains. 



'3a 



P / V P / P 

One can notice that I^a cancels the last term of Iia in Eq. ( |A17D . I^i, is harder to find. 
Substitution of Eq. ( [A 161 ) yields 



hb = -\ [s^diTiuG\' + iG'^d.Tu^ + Tius^AiG\' - G'^s'ATk 

= —-^s^{{di!Fik — 2TikAi)GQ^ + G}^{diJ^kj + '^^i^kj, 
1 

'2' 



kj 



-\s' {ViT,kG\' + C^^V.Tk, - {T,k. A} G\' - CT^ {Ai, T^,}) (A19) 
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We will split /3b further into two parts, I^b = hha + -^366, where 

I-iba = "2"^' {^l^ ikG^^ + G^qViT}:^ 

which, in the momentum representation, reads. 



i 

T4 



The rest of Isb is denoted as Jsab (which is the terms in Eq. ( |A19| ) that contains {J-", A}) 



and can be most conveniently written in the coordinate representation, 

hab = ^s' Ai} G',^ - Gi^ {Ai, J^kj}) = Is'^s' Ai} + {J^,k, Ai} s') G2 

Now when everything has been computed, one can collect all the terms. Many terms 
cancel each other (for example, I^ab cancels with Jit), and the equation for G has the form, 



k 



+ '^Vk:F,,Gs - ^p^p' {Vi^ikP^ - ViJ'.kp') Gs = Oig^) (A20) 

Notice that in the first term in the LHS of Eq. ( |A2(]| ), one can write p^V^G instead of p^V^G, 
since Dq is much smaller than D^. However, we will leave it as written in Eq. ( |A20| ) for 
future convenience. It can be check that the solution to this equation, up to the order of g^, 
can be written in the form. 



pV 'vIpI^IpI 2p2; ' p2 V- ^'^^ '''' J \2\^\d\v\ p2 

(A21) 

where we have introduced Gkin which satisfies the following equation. 
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- ip^V.G^,^ - P^^,u-Q^ = (A22) 



and made use of the fact that Gg depends on p via its module. One can recognize in Eq. 
( |A22|) the covariant form of the non-Abehan Vlasov equation. Making use of the exphcit 
form of Gsi Eq. (^) one can re-express G"^" in term of the new variable 5n'^{X), 

G^nX,p) = -ir''^5nl{X)2T,5{u^ - p2) (A23) 

where 8n^ satisfies the Vlasov equation in the non-covariant form, 

nil 

v^^D^5nl + g^¥.--^ = Q (A24) 

where n-p is the Bose-Einstein distribution function. Since E"" ~ g^T, one finds that 5n.p ~ 
g'^. Eq. (|A24D can be solved explicitly. The solution reads, 



cJu r 

(5np(t, x) = -g^r^ / du U''\t, x;t-u,x- uv)v'E\t - m, x - uv) (A25) 
Eqs. ( |A22| , |A23| , [7^^ ) completely determines the tree level G{x,y) up to the order of g^. 

2. Suppression of loop corrections 

In the previous section we have found tree- level G up to the order of g^ on any background 
field A (provided A ~ g^T). We have seen that some parts of the propagator can be 
attributed to the gauge invariance (the first term in Eq. ( A13|) ) and the deviation of the 



distribution function from thermal equilibrium {6n). Since small contributions, up to the 
order of g"^, are important for computing the damping term j in the stochastic equation, a 
question arises on the role of loop contributions to G: whether there is any loop diagram that 
contributes to the thermal propagator (and j) to the order of g^ or g'^. It is known that 
in certain cases the inclusion of loops leads to the Boltzmann collision term in the kinetic 
equation for 6n. So, the question we want to address is essentially whether the collision 
between hard particles change the dynamics of modes with g'^T spatial momentum. We will 
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present here a set of arguments showing that the corrections coming from the loops are are 
much smaller than the order of g we are interested in and can be neglected. 

In fact, according to the general formalism |13[, when one includes loops, Eq. ( |A11| ) will 
obtains a collision term in the RHS 

h + h + h = /col (A26) 

where 

Ico\{x, y) ~ I dz z)G{z, y) - G{x, y)) 



where S is the self-energy of the hard field. Since we are interested up to order in Eq. 
), we need to compute S up to this order, i.e. up to the three-loop level. 



To the one-loop level (Fig. |^), it is well-known that there is no contribution to Jcoi in 
the absence of background. This is because there is a hard particle cannot decay into two 
particles due to kinematic constraints. In the presence of the background, the analysis is 
more involved due to the complicated form of the propagator. However, essentially the same 
argument shows that there is no one- loop contribution to Jcoi. 

Two-loop graphs (for example, one in Fig. (0) where one of the internal propagator may 
be of G'^^ which corresponds to exchange of longitudinal gluon) are the lowest ones that 
generate the collision term in the absence of background. To find the same diagram on the 
background, one notice that the leading order of magnitude of the diagram is g"^, so one 
needs to compute these graphs to the leading and next-to-leading order. The propagator, 
up to the next-to-leading {g"^) order can be written as the that in the absence of background 
multiplied by the Wilson line, G{x,y) = U{x,y)Go{x,y), (we have shown that for spatial 
components of G, in fact it is also true for G^^). Now one can evaluate the graph in Fig. 
(§) by replacing each propagator line by UGq and the vertices by f^D. As the result S(x, y) 
on the background differs from the same quantity in the absence of the background by the 
Wilson line, 

T,{x,y) = U{x,y)i:{x,y)\A=o + 0{g^) 
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(the 0((y'^) corrections, roughly speaking, comes from the fact that U{x, z)U{z,y), where z 
is one of the internal vertices in Fig. |^ is not equal to U{x,y). However, the difference is 
of order g"^ which can be neglected). From that, one can show that the collision term Jcoi 
is the same as of on the absence of the background except for a Wilson line. However, it 
is known that the collision integral identically vanishes in thermal equilibrium (notice that 
Go is the propagator in thermal equilibrium). In particular, the two-loop contributions to 
the collision integral also vanishes in the absence of background, so the same is valid when 
background is present. 

Three- loop graphs are even simpler: they already have the marginal smallness g^, and 
to the leading order they can be computed in the absence of the background. Invoking the 
same argument about the vanishing collision integral in thermal equilibrium, we see that 
three-loop graphs have no contribution to the order of g^. 

Therefore, the propagator that has been found at tree level has no correction to the g^ 
order from loops. 

There is a simple (but not rigorous) way to understand the suppression of loop correc- 
tions. The collision integral Jcoi is proportional to 6n/T, where r is the relaxation time of 
the hard modes. Since 6n ~ g'^, and the relaxation time for hard modes is of order l/g'^T, 
the collision integral contains g^, which is much smaller than the order we are interested in 
(i.e. g^). 

APPENDIX B: PROPAGATORS AT LARGE SPATIAL SEPARATIONS 

To compute the noise kernel N^f^{x,y), one needs to find the propagator G{x,y) when x 
and y are separated by a distance of order {g'^T)~^. In contrast with the computation of the 
damping term j, here one needs to compute G only to the leading order. 

The propagator can be found most easily by making use of the Fourier expansion of the 
operator a\ Recall that in the absence of the background field (^4^=0), one can decompose 
the field operator a* into Fourier components as follows. 
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where a = 1, 2, e^(p) are the two transverse polarizations of a* that are perpendicular to p, 
and b"'" and W""^ are annihilation and creation operators for gauge bosons with polarization 
a and color a. 

In the presence of the background field A^, the equation for a*, up to contribution of 
order (we need to include, beside the leading term 0(1), also terms of order 0{g^), since 
the sub-leading, 0{g^) corrections may accumulate over the distance of order 0(((?^T)~^) to 
a large contribution), can be written as 

Dla'{x) = 

(cf. Eq. (^), note that F g'^). Eq. (^ is now modified to 

a^'{t, x) = I (^^^^ {f;\t, x)e-^IPl*+^PVJp)6jr + f^^t, x)e^lPl*-P-ej:(p)65r^) (B2) 

where /p*(t,x) satisfies the wave equation, 

D2(^/p(t,x)e-^IPl*+*P'') = (B3) 

To specify the boundary condition for /, one assumes that when t < the background field 
is turned off, A = 0, and demand that 

Let us check explicitly that, with the required accuracy, the solution to Eq. ( PB| ) is 

f;\t,^) = U^\t,^;0,^-vt) 

where U is the Wilson line. Indeed, we are interested in the leading and 0{g^) contributions 
to the equations, so 
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where the first term is the leading contribution, while the second term, containing one 
derivative of a slowly- varying U, is of order 0{g'^). Taking one more derivative, one finds, 

where we have neglected the contributions like with two derivatives of U since it is of order 
0((7'^). Now, as we have = and p^D^U = 0, / satisfies the wave equation Eq. ( |B^j| ) with 
the accuracy that we are interested in, i.e. up to 0{g^) corrections inclusively. 

Now making use of the fact that, in thermal equilibrium, (6p&p') = (27r)^2|p|r2p5(p — p'), 
where rip is the Bose-Einstein distribution function, the propagator can be written in the 
form, 

GZi^^y) = {U^)f,\y)T {^'^ - ^) (npe-ipi(*-*')+-(--) + 

+ (l + np)e'|p|(*-*')-'P("-y)) (B4) 



APPENDIX C: ONE-DIMENSIONAL LANGEVIN EQUATION WITH NOISE 
AND DAMPING DEPENDENT OF COORDINATE: ASYMMETRIC 

DISCRETIZATIONS 

In this appendix we will consider the following one- dimensional stochastic equation 

^{x)x + V\x)=i (CI) 

(e(t)e(0)=2T7(x)5(t-t') (C2) 

Let us first neglect the V' term in Eq. ( |C1|) . Let us first choose the simplest discretization 
of Eqs. ( |CTlJC2|) . Consider the following finite-difference equation. 



. ^^ 2Tj(x(t)) 
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Suppose at time moment t the particle's coordinate is y, then the probabihty of finding, at 
time moment t + At, the particle at coordinate x is 



Pit, y — > t + At, x) = \ — — — exp — 7^77— [v — x) 
V ' y ^ V 47rTAt ^ V 4TAt ' 

Thus, the evolution of p{t, x) is described by the following equation. 



pit + At,x) = JdyPit,y)j£^^exp(-^^{y-xr\ (C3) 



Now one notes that when At is small, the integration in Eq. ( |U3D is saturated mostly by the 
region of small y — x. Therefore, one can expand the integrand in Eq. (|C^ ) around y = x 
and take the integral. One finds, that in the limit of At 0, the evolution of p(t, x) follows 
the following partial differential equation (the Fokker-Planck equation). 

When the V term in Eq. ( |CT| ) is taken into account, the Fokker-Planck becomes. 

Now let us choose another discretization of the Langevin equation. Instead of taking the 
value of 7 at x(t), let us take its value at x{t + At). In other words, consider 

7(x« + A<))^(^±M^ = ^W (C5) 
«W«* + „A()> = ?lM^*„o (C6) 



Eqs. ( |C5|JC6|) may present some difficulty for their understanding: to generate the noise 



at the time moment t, according to Eq. (|Cq), one needs to know the coordinate at the 
time moment t + At. We will understand Eqs. (|C5| , |C6|) as we did for Eqs. (|1^,|^), i.e. by 
considering a statistical ensemble of all possible path x(t) where each path is associated with 
a weight. 



p[x{t)] ~ exp I - X] 



4T7(nAt) 



7(nAt) 



x(nAt) - x{{n - l)At) 



1 2> 



At 
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In this understanding, the probabihty of going from the point y at t to the point x at time 
t + At is 

P{t, y-^t + At,x) = C exp (-^(2/ - xf^ (C7) 

where C is the normahzation constant defined by the condition of probabihty conservation, 

JdxP{t,y ^ t + At,x) = 1 

The Fokker-Planck equation can be derived similarly with the previous method of discretiza- 
tion, now the result reads, 

— -T—(——h^p)]+ — (—p] (C8) 
dt dx dx ^ ' j (9x \ 7 / 

where the contribution from the term V in Eq. ( |Cl|) has been also taken into account. 
It is clear that this equation is different from Eq. ( P4|) , derived using another scheme of 
discretization. Therefore, we see that the motion of the particle is not uniquely defined by 
Eqs. (|CT],0). 

APPENDIX D: ONE-DIMENSIONAL LANGEVIN EQUATION WITH NOISE 
AND DAMPING DEPENDENT OF COORDINATE: A SYMMETRIC 

DISCRETIZATION 



In this Appendix we will consider the same stochastic equation as in Appendix Eqs. 
( Pl| , |C^ ) . In Appendix |C| we have seen that the Fokker-Planck equation depends on the the 



way one "regularizes" the stochastic dynamics. Here we will use the following prescription. 
We add a small second time derivative term to the LHS of Eq. 

ex + '-fx + V'{x) = ^(t) 



and understand this equation and Eq. ( |U2D as the limit A — of the following finite- 
difference equations. 
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^.(i + At) - 2g) + .(t - At) ^ ^^^^^^^.jt + A(,^(t - At) ^ ^^.^^^^^^ ^ ^^^^ 

2T 

(e(t)e(i + nAi)) = — 7(x(t))5„o (D2) 

Notice that x is approximated with the accuracy of (At)^, in contrast with the prescriptions 
in Appendix |C| where the accuracy is only At. The hmit A ^ should be taken before the 
limit e —* 0. In other words, we assume that A ^ e, while both are small. It is useful to 
introduce the momentum, 

_ x{t) - x{t - At) 
~ At 

Eq. ( pi| ) now can be recasted into the form. 

Neglecting small corrections, the evolution in the {x,p) space is described by the following 
couple of equations, 

x{t + At) = x{t) + p{t + At) At 

p{t + At) =(l- ^^^^At\ p{t) - ^V'{x{t)) + ^at) (D3) 



These equations, together with Eq. ( |D2| ) completely defines the evolution of the stochastic 
system. The state of the system at any time slice t is described by a couple of variable 
{x,p). For a given state, one generates a Gaussian noise ^ having the mean square equal to 
2T'j{x{t))/At and run the system one time step forward using Eqs. ( p3| ) . 

Our strategy is as follows. Denoting by p(t,x,p) the probability distribution of the 
particle in the {x,p) at time moment t, we will find the equation describing the evolution 
of p{t,x,p) with time. After that, we will see that at time scales large compared to e 
this equation reduced to the equation for the distribution function in the x-space alone, 
p{t,x) = Jdpp{t,x,p). 

Assuming ^ is Gaussian, from Eqs. (|D3|) one finds, 
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p{t + At,x,p) = J dydqp{t,y,q)5{x - y - pAt) 



X exp 



p-il- -At)q + —V'{x 



J dy dq p{t, x - pAt, q) ^^Jj,^^^ exp 



4T7At V V e / ^ e 



In the limit At ^ 0, the integral is saturated by the integration region where q is closed to 
p, where one can replace by its Taylor expansion around p, p(g) = p(p) + (g — p)p'{p) + 
— pYp'^p) + ■ ■ ■• In this limit, one finds that p{t,x,p) satisfies the following diffusion 
equation, 

dt dp'^ ^ dx^ e ~^ ^ dp^ ^ ^ 

In all our discussions above we have kept e finite. Now let us take the limit e and try 
to find the evolution of the distribution function on time scale of order 1. Assuming that 
the p is of order e~^/^ (we will check this assumption a posteriori), if one keeps only leading 
on e terms in the RHS of Eq. (|D4|) , the latter reduces to 

dp T7 d"^ p 1 dp 7 
dt dp"^ dp 

The equation has solutions in the form. 



p{t,x,p) = ^ 



27rr 

where a is an arbitrary function of x. For this solution the typical value of p is obviously of 
order e~^/^. Now, let us look for the solution to Eq. ( [D4| ) in form of the following series, 

p{t,x,p) = po{t,x,p) + pi{t,x,p) + p2{t,x,p) H (D5) 

where the expansion parameter is e^/^, pi/po ~ P2I Pi ~ e^''^, and 



50 



where a(t,x) is a function that varies on the times scale of order e°, d/a ~ 1. Substituting 
Eq. (p5| ) to Eq. ( p4| ) and putting together terms of the same order on e, one finds the 
following equations for a, pi and p2, 

1 (Z^ +p^+p^ + (Yl^- p^] = (D7) 
eye dp"^ dp J \ e dp dx J 

-epy2Tda 1 (Td'^p2 dp2 \ (V'dpi dpi\ 

Substituting po in Eq. (|D6|) to Eq. ( P7|) , one obtains the following equation, 

eye dp"^ ^ dp ^ y ^ \dx T j 

from which one can find pi, 



Substituting this expression for pi to Eq. ( p8| ), one find p^ 

and the equation that a must satisfy, 

— - T— - f — + — 

According to Eqs. (P5|) and (P6|) , a can be considered as the distribution function of particle 
with respect to the coordinate alone. In other words, with the accuracy e^/^ 

a{t,x) = jdpp(t,x^p) 

Re-denoting a(t, x) as p{t,x), we find, finally, the Fokker-Planck equation, 

—p(t x)- — i^—p{t x) \ + — { ^'^^\ {t x)\ (D9) 
dt ' dx \7(x) dx ' y dx \ 7(x) ' y 
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FIGURES 
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FIG. 1. The Schwinger-Keldysh Close-Time-Path contour 




FIG. 2. A graph with one hard loop 




FIG. 3. A graph with two hard loops 
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